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A notorious difficulty in the covariant dynamics of classical charged particles subject to non-local 
electromagnetic (EM) interactions arising in the EM radiation-reaction (RR) phenomena is due to 
the definition of the related non-local Lagrangian and Hamiltonian systems. As a basic consequence, 
the lack of a standard Lagrangian/Hamiltonian formulation in the customary asymptotic approxi- 
mation for the RR equation may inhibit the construction of consistent kinetic and fluid theories. In 
this paper the issue is investigated in the framework of Special Relativity. It is shown that, for finite- 
size spherically-symmetric classical charged particles, non-perturbative Lagrangian and Hamiltonian 
formulations in standard form can be obtained, which describe particle dynamics in the presence 
of the exact EM RR self-force. As a remarkable consequence, based on axiomatic formulation of 
classical statistical mechanics, the covariant kinetic theory for systems of charged particles subject 
to the EM RR self-force is formulated in Hamiltonian form. A fundamental feature is that the non- 
local effects enter the kinetic equation only through the retarded particle 4-position. This permits, 
in turn, the construction of the related fluid equations, in which the non-local contributions carried 
by the RR effects are explicitly displayed. In particular, it is shown that the moment equations 
obtained in this way do not contain higher-order moments, allowing as a consequence the adoption 
of standard closure conditions. A remarkable aspect of the theory is related to the short delay-time 
asymptotic expansions. Here it is shown that two possible expansions are permitted. Both can be 
implemented for the single-particle dynamics as well as for the corresponding kinetic and fluid treat- 
ments. In the last case, they are performed a posteriori, namely on the relevant moment equations 
obtained after integration of the kinetic equation over the velocity space. Comparisons with the 
literature are pointed out. 



I. INTRODUCTION 



An open problem in relativistic theories is related to the Hamiltonian description of particle dynamics for which 
non-local interactions typically occur. In this regard, a basic difficulty which is usually met is the lack of a Hamiltonian 
formalism for non-local Lagrangian systems. In fact, for arbitrary non-local Lagrangians it is generally impossible to 
define the notion of Legendre transformation Q. As a consequence even the phase-space itself may not be well-defined. 

Most approaches to the construction of a Hamiltonian formalism for non-local first-order Lagrangians have tried 
to change the functional part of the Euler-Lagrange equations In principle this delivers infinite-order Euler- 

Lagrange equations and a corresponding infinite-dimensional phase-space. As an alternative, a finite dimensional 
phase-space can be recovered by introducing appropriate asymptotic approximations, i.e., truncating the expansion 
of the Lagrangian in terms of finite-order derivatives [3, @ . 

A typical situation of this kind occurs for the relativistic equation of motion for single isolated charged particles, 
subject both to external and self EM forces, namely the radiation-reaction (RR) equation. There is an extensive 
literature devoted to this subject, most of which dealing with point charges. As remarked by Dorigo et al. 0], 
customary formulations based either on the LAD [8hTo| or LL [ill ] equations are asymptotic, i.e., obtained by means of 
asymptotic expansions of different sort. In particular, as a consequence it follows that the LAD equation is represented 
by a third-order ODE, so that it does not admit a Hamiltonian formulation in the customary sense [HI, [l3j]. The 
LL, instead, is intrinsically non-variational, although it is a second-order differential equation, being obtained by 
means of a one-step "reduction process" from the LAD equation p}. As a consequence, the LAD equation does 
not define a dynamical system in the customary sense, since it requires, for non-rotating particles, a 12-dimensional 
phase-space involving also the particle acceleration. Therefore, for different reasons, both the LAD and LL equations 
are manifestly non-Hamiltonian. In particular, for the LL equation, this implies that the corresponding phase-space 
volume is not conserved. Moreover, within these treatments particles are treated as point-like, so that non-local EM 
effects produced by the RR self-interaction may remain undetermined. 

Fundamental problems arise when attempting to formulate classical statistical mechanics (CSM) for systems of 
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relativistic charged particles based on the LAD or LL equations. In fact even the proper axiomatic formulation of the 
rclativistic CSM for radiating particles is missing. This requires the precise identification of the corresponding phase- 
space and the definition of an invariant probability measure on this set. For a system of charged particles subject 
solely to an external EM field and the RR self- force this involves the construction of a Vlasov kinetic treatment. In 
this regard, important issues concern: 

1) The lack of a standard Hamiltonian formulation of relativistic CSM based on such asymptotic equations, which 
implies the lack of a flow-preserving measure. This feature is shared by both the LAD and LL equations. 

2) The proper definition of a phase-space. The problem is relevant for the LAD equation. In fact, although the 
construction of kinetic theory is still formally possible [3, HH j the corresponding fluid statistical description seems 
inhibited. 

3) The explicit dependence of the kinetic distribution function (KDF) in terms of the retarded EM self 4-potential are 
excluded. In fact, in the LAD and LL approximations the self-potential does not appear explicitly (see for example 
Refs.flU [l7|). Indeed, within the point-particle model, underlying both treatments, the retarded self-potential is 
divergent. 

On the other hand, for the fluid treatment: 

1) The precise form of the fluid closure conditions may depend on the approximations adopted in the kinetic 
description for the representation of the EM RR self- force. An example-case is provided by Ref.[17| where relativistic 
fluid equations are obtained based on the LL equation. As a result, it was found that, with the exception of the 
continuity equation, all moment equations involve higher-order fluid moments associated to the RR self-force. It is 
unclear whether this is an intrinsic physical feature of the RR interaction or simply a result of the approximations 
involved. 

2) The fluid fields may in principle depend implicitly on the EM self 4-potential. In the framework of the LL equation 
it is unclear how such an effect can be dealt with. However, the treatment of such effects seems to present objective 
difficulties. In fact, in principle non-local effects might arise in this way in which retarded velocity contributions appear 
in the kinetic equation. In such a case the explicit construction of fluid equations would be ambiguous (and might 
involve an infinite set of higher-order moments) . 

The interesting question is whether these difficulties can be overcome in physically realizable situations, namely for 
exactly solvable classical systems (of particles) for which the relativistic equations of motion are both variational and 
non- asymptotic. The prerequisite is provided by the possibility of constructing an exact representation for the RR 
equation for a suitable type of classical charged particles. In the past, their precise identification with physically- 
realizable systems has remained elusive because of the difficulty of the problem. However, as recently pointed out 
(see Tessarotto et al. [l8| and Cremaschini et al. [ijjj . hereafter referred to as Paper I) in the framework of special 
relativity an exact variational RR equation can be obtained for classical finite-size charged particles. This refers to 
particles having a finite-size mass and charge distributions which are quasi-rigid, non-rotating, spherically symmetric 
and radially localized on a spherical surface 9f2 having a finite radius a > (see [2(| and the related discussion in 
Paper I). In this formulation, contrary to the point-particle case, the retarded EM self 4-potential is well-defined, 
namely, it does not diverge, and can be determined analytically. As shown in Paper I, it follows that the RR 
equation is variational and the corresponding Hamilton variational functional is symmetric with respect to the non- 
local contributions. The latter are due to the retarded EM self interaction arising from the finite spatial extension of 
the charge distribution. As a consequence, the resulting exact RR equation is a second-order delay-type ODE which 
admits a Lagrangian formulation in standard form (see discussion below). Furthermore, under suitable conditions, 
the same equation defines a classical dynamical system (RR dynamical system). 

In this paper we intend to prove that, based on the results of Paper I, the RR dynamical system admits also a 
Hamiltonian representation in terms of an effective non-local Hamiltonian function H e *f . This implies that the exact 
RR equation can also be cast in the equivalent standard Hamiltonian form represented by first-order delay-type ODEs 



with y = (r M , P^) denoting, in superabundant variables, the particle canonical state which spans the eight-dimensional 
phase-space T = r r x T,, where T r and T u are respectively the Minkowski M 4 -configuration space and the 4- 
dimensional velocity-space, both with metric 7/ M „ = diag (1, — 1, — 1, — 1). Remarkably, here it is found that the 
Hamiltonian structure can be retained also after the introduction of a suitable short delay-time approximation of the 
RR force. The result is an intrinsic feature of the extended particle model adopted in the present treatment. 

As a consequence, the statistical description of the RR dynamical system follows in a standard way. In particular, 
here we report both the exact and asymptotic kinetic and fluid formulations. These are developed for collisionless 



dr^ 
ds 

ds 



Q H eff 



(1) 



Qfieff 

dr^ 
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relativistic plasmas in the Vlasov-Maxwell description, including consistently the contribution carried by the RR 
self-field. Applications of the theory here developed concern: 

1) The kinetic and fluid treatments of relativistic astrophysical plasmas observed, for example, in accretion disks, 
relativistic jets, active galactic nuclei and mass inflows around compact objects. 

2) The kinetic and fluid treatments of laboratory plasmas subject to ultra-intense and pulsed-laser sources. 

A. Goals and scheme of the presentation 

In detail, the plan of the paper is as follows. In Section 2 we introduce the basic definitions concerning the 
Lagrangian formulation for non-local interactions and the concept of effective non-local Lagrangian function and the 
related covariance properties (THM.l and Corollary). In Section 3 we provide an analogous generalization which 
permits to introduce the notion of non-local Hamiltonian formulation (THM.2 and Corollary). In particular, within 
the framework of the theory here developed, the standard form for the local Legendre transformation is retained, while 
the concepts of effective canonical momenta and effective Hamiltonian function are introduced. Then, in Section 4 the 
general case of a rotating finite-size and spherically-symmetric charged particle is discussed as a physical example of 
non-local interaction. It is shown that the corresponding Hamilton variational functional satisfies the requirements of 
THMs.l and 2 (see THM.3) and therefore admits both Lagrangian and Hamiltonian formulations in standard form. In 
Section 5 the theory is applied to the specific case of a non- rotating particle. As a result, based on the analytical results 
of Paper I, the corresponding variational and effective Hamiltonian formulations are presented (THM.4). This provides 
the explicit form of the effective Hamiltonian function and a parameter-free representation for the retarded EM RR 
self- force. Then, in Section 6 a Hamiltonian asymptotic approximation of the RR equation is developed (THM.5), 
based on the retarded-time expansion holding in the short delay-time ordering. The approximation overcomes basic 
inconsistencies of the LAD and LL treatments, applying in the case of point-particles. As a consequence, in Section 
7 the relativistic kinetic theory for a collisionless plasma with the inclusion of the EM RR effect is formulated in 
Hamiltonian form. The use of superabundant canonical variables allows a precise identification of the phase-space 
and the consequent axiomatic formulation of CSM with non-local EM RR interactions. This is based on the notion of 
invariant probability measure in such a setting. As a result, a relativistic Liouville-Vlasov kinetic equation is proved to 
hold for the KDF (THM.6). This permits to achieve a Vlasov-Maxwell description applicable to relativistic plasmas, 
in which the RR interaction is consistently taken into account also in the Maxwell equations both for the external 
and self EM fields. In Section 8 the corresponding fluid fields and fluid moment equations are computed in terms 
of 4-velocity integrals, which retain the standard conservative Eulerian form as in the absence of RR effects. The 
existence of both explicit and implicit non-local contributions arising from the RR effect in the fluid equations is 
discussed. It is shown that the former are associated to the EM force acting on the fluid, while the latter enter in 
the definition of the fluid fields through the effective momenta. In particular, the explicit dependence of the KDF on 
the retarded EM self 4-potential is discussed and an asymptotic estimation of the implicit contributions is presented. 
In Section 9 a Lagrangian formulation of the fluid equations is derived, which allows one to introduce an explicit 
parametrization of the non-local RR terms carried by the EM self 4-potential and the EM RR force. It is shown that 
the exact fluid equations with the inclusion of the RR interaction are of delay-type. Section 10 deals instead with 
the development of asymptotic approximations of the moment equations. This is motivated by the requirement of 
reducing the exact non-local fluid equations to a local form. Different asymptotic approximations are obtained for the 
non-local terms of the RR effect, based both on short-delay time expansions (THM.7) and an iterative procedure which 
holds under the assumption of weak RR self- force (in comparison with external EM and pressure forces). A detailed 
analysis of the basic physical properties of the kinetic and fluid treatments obtained here and a comparison with 
previous literature is reported in Section 11. Concluding remarks are presented in Section 12. Finally, in Appendix 
A a Green-function approach is developed for the calculation of the EM self 4-potential, while in Appendix B the 
connection with non-canonical representations is provided for the relativistic kinetic theory. 

II. NON-LOCAL LAGRANGIAN FORMULATION 

The natural mathematical apparatus for an abstract description of Lagrangian and Hamiltonian mechanics is that 
of variational principles, whose methods have been studied for a long time by mathematicians and can be found in 
the textbooks. Nevertheless, actual problems of interest in classical relativistic dynamics involving the treatment 
of non-local interactions have escaped a solution. In particular, in the literature the prevailing view is that, while a 
non-local variational formulation is possible, a corresponding Hamiltonian representation is generally excluded. In the 
following we intend to point out that for a particular class of non-local Lagrangian systems the problem can be given 
a complete solution. The latter correspond to variational problems in which the variational functional is symmetric. 
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To this end, in this section we briefly recall basic notions holding for local and non-local Lagrangian systems. This 
task represents a necessary prerequisite for the establishment of a corresponding Hamiltonian formulation and for the 
subsequent investigation of the Hamiltonian dynamics of finite-size charged particles with the inclusion of the RR 
self-force. 



Definition #1 - Local and non-local Lagrangian systems. 

A local (respectively, non-local) Lagrangian system is defined by the set {x, L} such that the following conditions 
are satisfied. 



1. x = (r^ (s) , dr d ^ = r \ is the Lagrangian state spanning the Lagrangian phase space C 

2. The Lagrangian action functional S is a 4-scalar of the form 



v2N 



S2 

dsL, (2) 

si 



with L to be referred to as variational Lagrangian function. In particular, the functional dependencies of S and 
L are respectively of the form: 



• S = So (r) and L = Lq (r, ^) for local systems; 

• S = Si (r, [r]) and L = L\ (r, 4^, [r]) for non-local systems, with [r] denoting non-local dependencies. 

3. In the functional class 

{r"} = K(s) : r"(«<) = rf, for i= 1,2, si,s 2 € /, 

with si < s 2 and r^(s) G C 2 (7)} , (3) 

the synchronous variations (5r M (s) are considered independent and vanish at the endpoints r^(sj) = . Hereafter 
6 denotes, as usual, the Frechet functional derivative. For a synchronous variational principle the interval ds is 
such that Sds — and is subject to the constraint 

ds 2 = g^dr" (s) dr v (s) , (4) 

where r M (s) are the extremal curves. 

4. The Lagrangian action ([2J admits a unique extremal curve r A '(s) such that, for all synchronous variations <5r M (s) 
in the functional class ([3]) the Hamilton variational principle 

SS = (5) 

holds identically. For non-local systems the non-local Lagrangian must be suitably constructed in such a way 
that the extrema curves r M (s) satisfy the constraint (j4j. 

In particular, for local systems the extremal curves of Sq are provided by the Euler-Lagrange (E-L) equations 

H = F,(r)L = 0, (6) 

where, for an arbitrary set of Lagrange coordinates g M , F^q) denotes the E-L differential operator 

T-, i \ d d d ,_. 



and q = fa^ (s) . 



On the other hand, for non-local systems the extremal curves of the functional 5i are provided by the Euler- 
Lagrange equations 



SSi = SSi 
8r^ Sr^ 



5Si 



0, (8) 



where |%| r . and 4% 



carry respectively the contributions due to the local and non-local dependencies. 
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Definition #2 - Non-local Lagrangian systems in standard form. 

A non-local Lagrangian system {x, Li} will be said to admit a standard form if the variational derivative (JSJ yields 
the E-L equations in the standard form: 

§ + ^-W e// =0, (9) 

with 

L eff = L ef f(r,^,[r]) (10) 



denoting a suitable effective non-local Lagrangian. 



On the base of these definitions, the following theorem holds. 



THM.l - Non-local and Effective Lagrangian functions 

Given validity of the definitions #1 and #2, it follows that: 

Tli) The non-local Lagrangian L\ and the effective Lagrangian L e ff are generally different, namely 

Li^L eff . (11) 

Tl 2 ) If Si (r, [r] ) admits the general decomposition 

S 1 (r,[r]) = S a (r)+S b (r,[r]), (12) 

with S a {r) = f^ 2 dsL a (r, ^) and Sb(r, [r]) = J^ 2 dsLf, (r, [r]), and moreover Sb(r, [r]) defines a symmetric func- 
tional such that 



S b (r,[r])^S b ([r},r), (13) 

then the effective Lagrangian L e ff is related to the variational non-local Lagrangian L\ = L a + L b as 

L ef f = L a + 2L b = Li + L b . (14) 

Proof - Tli) The proof is an immediate consequence of Eqs.® and (|9]). In fact, by definition the E-L differential 
operator F^ir) is a local differential operator that is required to preserve its form also for non-local systems. On the 
other hand, the variational derivative © is different from Hence, in order to write the E-L equations associated 
to the non-local function L\ in standard form, a suitable effective Lagrangian £ e // must be introduced, which must 
differ from L\ and be expressed in such a way that the non-local dependencies contained in L\ can be equivalently 
treated by means of F^{r). 

TI2) The proof follows by inspecting the general definition ([S]). In this case, in view of the symmetry property (|13l) . 
it follows manifestly that 

5Sx _ SSi 

Then, by comparing this relation with the definitions both of the E-L differential operator (JTJ) and the standard form 
representation of the E-L equations ©, it follows that the effective Lagrangian L e ff takes necessarily the form given 
in Eq. (fl"4]) . This completes the proof of the statement. 
Q.E.D. 



SS 1 








Sr^ 





(15) 



A basic consequence of Definition #2 and THM.l concerns the covariance property of the E-L equations ©. The 
result is stated in the following Corollary. 



Corollary 1 to THM.l - Covariance of the E-L equations for arbitrary point transformations. 

The Euler- Lagrange equations 0) are covariant with respect to arbitrary point transformations 



r^ £^(r) 



(16) 
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represented by a diffeomeorphism of class C k , with k > 2, which requires they are of the form 

F MLeff = ^F u (r)L eff = 0, (17) 

with L e ff denoting 

1*/ (-.&«)■ (i« 

j4s a consequence, Eq. satisfies also the covariance property with respect to arbitrary infinitesimal Lorentz trans- 
formations (Manifest Lorentz Covariance). 

Proof - The Euler-Lagrange equations © are by definition covariant provided the variational Lagrangian 
L\ (r, j-)H) is a 4-scalar (as it is by construction). Then, it is sufficient to represent the Lagrangian action in 
terms of the Lagrangian coordinates q^, yielding 

&(q,[q])=Si(rM, (19) 

with Si(q, [q]) denoting the transformed action 

& 1 (q i [q])=J''dsT 1 (q,^,[q]) (20) 

and L\ denoting the transformed variational non-local Lagrangian. Hence, the Hamilton variational principle 
SSi(q,[q]) = yields precisely the E-L equations (TTT)) . This proves the statement. The covariance property of 
Eqs.Q with respect to point transformations (|1 6[) includes, as particular case, Lorentz transformations. Therefore, 
Eqs.© are also Manifestly Lorentz Covariant (MLC). 
Q.E.D. 

We notice the following remarkable features of this treatment: 

1) In general, in absence of any kind of symmetry, a non-local Lagrangian system does not admit a standard form 
representation in terms of the effective Lagrangian L e j f [l[ . 

2) As shown in TI2, the possibility of getting an explicit relationship between L\ and L e ff is a consequence solely 
of the symmetry property (II 3[) of the functional Sb- This also proves the existence of £ e // and, as a consequence, of 
the standard form representation for non-local systems satisfying Eq. (fT3|) . 

3) The symmetry assumption (|13[) can be effectively realized in physical systems. As it will be shown below, this 
condition is satisfied by the variational functional which describes the dynamics of finite-size classical charged particles 
with the inclusion of the RR effects associated to the interaction with the EM self-field. 



III. NON-LOCAL HAMILTONIAN FORMULATION 

In this section we deal with the basic features concerning the Hamiltonian formulation for non-local systems which 
admit a variational treatment in terms of non-local Lagrangian functions. This requires the introduction of the 
following preliminary definitions. 

Definition #3 - Local and non-local Hamiltonian systems. 

A local (respectively, non-local) Lagrangian system {x, L} is said to admit a local (non-local) Hamiltonian system 
{y = (r^jPf^jH} provide the following conditions are satisfied. 

1. The variational Hamiltonian H is defined as the Legendre transformation of the local (non-local) variational 
Lagrangian L 

dr^ 

H = P/i — — L, (21) 



with 

dL 
d 



V, = t^tt (22) 



being the corresponding canonical momentum, with corresponding action functional 5*^ = ds [p^ — if] . 
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2. It is assumed that H is respectively of the form: 



• H = Hq (r,p) for local systems; 

• H = (r, p, [r]) for non-local systems, namely it is a local function of (r,p) and a functional of [r]. 



The corresponding Hamilton action Junctionals are denoted respectively as 



SH (r 7 p) = ds 

J ax 



for local systems, and as 



for non-local systems. 
3. In the functional class 



ShA^P, M) = / ds 



dr^ 
Pv-nr-Hi 



(23) 



(24) 



{y = (^ 5 P/i)} = {y( s ) ; y{ s i) = y*, for * = M, si,s 2 e 7, 

with si < s 2 and y(s) S C 2 (/)} 



(25) 



the synchronous variations (6r fl (s),Sp IJl (s)) are all considered independent and vanish at the endpoints y(«i) 
yi. By assumption, synchronous variations imply that <Ms = 0, with the interval ds satisfying the constraint 



ds 2 = g^dr" (s) dr v (s) . 



(26) 



where (s) are the extremal curves. 
4. The modified Hamilton variational principle 



ss 



H 







(27) 



with variations (5r^(s), 5p^(s)) is equivalent to the Hamilton principle (|5|), i.e., it yields the same extremal 
curves in the functional class {y}. 

In particular, for local systems the extremal curves of Sh can be cast in the standard Hamiltonian form as 
first-order ODEs 



SS Ho dr^ dH 



bPv 



ds 



9pn 



\r»H 



o 



(28) 



SS 



Ho 



dp^ _ dH 



[Pn,H ] , 



(29) 



8r^ ds dr^ 

where the customary Poisson bracket formalism has been used. 

On the other hand, for non-local systems the extremal curves of the functional Shx are provided by the set of 
first-order ODEs 

5S Hl 



SPl 



0. 



(30) 



where 



and 



SS H 
8[r^ 



ss Hl _ ss Hl 

Sr^ Sr^ 



6S 



Hi 



S 



(31) 



carry respectively the contributions due to the local and non-local dependencies. 
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Definition #4 - Non-local Hamiltonian systems in standard form. 

A non-local Hamiltonian system {y, H±} will be said to admit a standard form if the extremal first-order ODEs 
(|30[) and (|31| can be cast in the standard Hamiltonian form in terms of the effective canonical momentum P^ and 
Hamiltonian function H e ^ as 



SS Hl = dr^ = dHeff 



ds dP„ 



. H eff] 



ss 



SS 



Hi 



5S 



Hi 



5 frH 



dP» _ dH eff 
dr^ 



ds 



(32) 
(33) 



Here both H e f f — H e f f (r, P, [r] ) and P M must be dehned in terms of the effective Lagrangian function introduced in 
Eq. (|9j) respectively as 



H -P dr " T 



and 



Pa = 



OLeff 



(34) 



(35) 



From this definition it follows that, if the non-local Hamiltonian system {y,Hi} admits a standard form, then the 
Poisson bracket representation holds for He// and P u . 

The following theorem can be stated concerning the relationship between Hi and H e f f . 

THM.2 - Non-local and Effective Hamiltonian functions 

Given validity of the definitions #3 and #4 an d the results of THM.l, if Sh x (r,p, [r]) admits the general decom- 
position 



Sh! (r,P, H) = S Ha (r,p) + S Hb {r,p, [r]), 



with 



S H a (r,p) 



S Hb {r,P, M) = 



ds 



ds 



dr* 1 



dr^ 



where the canonical momenta p a ^ L and are defined respectively as 

dLa 



Pa,, 



Pbfj, 



ds 



d 



and moreover SH b (r,p, [r]) defines a symmetric functional such that 

SH b {r,p, [r]) = S Hb {[r] ,p,r), 
then the effective Hamiltonian H e ff is related to the variational non-local Hamiltonian Hi = H a + Hi, as 

Heff = 

where, by definition 



+ 2H b = 




dr" 






Li, 


dr» 




Pati lls~~ 


- L ai 


dr^ 






-L b . 



(36) 

(37) 
(38) 

(39) 
(40) 

(41) 
(42) 

(43) 
(44) 
(45) 
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Proof - The proof follows from THM.l and by invoking the general definitions (|32|) and ([33)) . In fact, in view of 
the symmetry property (flTj) , it follows manifestly that 

(46) 

[r] 

Then, by comparing this relation with the definitions ([34 |) - ([35l) for the standard Hamiltonian form and using Eqs. ([43|) - 
(|45p . from the analogous results in THM.l which concerns the relationship between L\ and L e tf in the symmetric 
case, Eq. (|42| is readily obtained. 



SSh! 


_ ss Hl 


, 6S Hl 


._ ss Ha 


+ 2 — 


Srt* 


Srt* 


[r] S M 







Q.E.D. 

Finally, as a basic consequence of Definition #4 and THM.2, the following Corollary can be stated concerning the 
covariance property of the Hamilton equations in standard form. 

Corollary 1 to THM.2 - Covariance of the Hamilton equations for arbitrary point transformations. 

The Hamilton equations H3S^) - [33\) in standard form are covariant with respect to arbitrary point transformations 

r M _> 9 M( r ) (47) 
represented by a diffeomeorphism of class C k with k > 2, which requires they are of the form 



dq» 
ds 



dH, 



dP, 



dP {q) ^ 
dH eff 



ds 



dq^ 



I^Heff 



P(q)f^H e ff 



(48) 
(49) 



with H e ff denoting 



Heff (q,P,[q]) = H eff (r,P,[r]) 



(50) 



and P( q )^ being the transformed canonical momentum. As a consequence, Eqs. |^ff| ) and |^ff| ) satisfy also the covariance 
property with respect to arbitrary infinitesimal Lorentz transformations (Manifest Lorentz Covariance). 

Proof - In fact, for an arbitrary point transformation of the type (|47[) . the corresponding transformation for the 
momenta P v is 



which yields 



Hence, it follows that 



dP, 



dPa 



dP, 



dP v 



dq» 
~dr»' 
dr v 
dqt' 



(51) 

(52) 
(53) 



which implies 



dq» 
ds 



dP, 



ds 



dH, 



eff 



dP iq)fi 

dH e ff 



dq» dr v 
dr v ds 



t dHeff 



dr" dP v ' 
dP( q)li dP v _ Or" OHeff 



dq» 


dP v 


ds 


dHeff 


dq» 


dHeff 


dP(q)» 


dr v 


dP v 


dHeff 


dr v 


dHeff 


dq^ 


dqf 1 


dr v 



(54) 
(55) 

(56) 
(57) 
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where H e ff is defined in Eq. ([50| above. Therefore, the Hamilton equations in standard form for the Lagrangian 
coordinates and the canonical momenta -P( g ) M are respectively covariant [Eq. (j56")l ] and controvariant [Eq. ([57)l ] with 
respect to the point transformation (|47|) . This is true also for arbitrary infinitesimal Lorentz transformations, which 
proves the MLC of Hamilton equations Eqs. (|4"5)l and f4"9")l in standard form. 
Q.E.D. 



IV. AN EXAMPLE OF NON-LOCAL INTERACTION: THE CLASSICAL EM RR PROBLEM 

A crucial issue of the present investigation concerns the possible existence of physical systems subject to non-local 
interactions whose dynamics can be consistently described in terms of a variational action integral and which admit 
at the same time both Lagrangian and Hamiltonian formulations in standard form. In this section we prove that the 
EM RR problem for classical finite-size charged particles represents a physical example of non-local interactions of 
this kind. The reason behind the choice of considering extended particles is the necessity of avoiding the intrinsic 
divergences of the RR effect characteristic of the point-particle model. 

In fact, consider the general form of the Hamilton action functional for the variational treatment of the dynamics 
of an extended charged particle in presence of an external EM field and with the inclusion of the RR self-interaction. 
This can be conveniently expressed as follows: 

S l (z, [z]) = S M (z) + S ( c ext) (z) + S ( c self \ Zl [z]), (58) 

where Sm, S^^and Sq b1 ^ are respectively the contributions from the inertial mass and the EM coupling with the 
external and the self fields. In particular, denoting by j( seZ /)^(r) the particle 4-current density generated by the 
particle itself and observed at a 4-position r, the two coupling action integrals are provided by the following 4-scalars: 



1 f 2 

S ( c xt) {z) = — J d 4 rA( ext ^(r)jl sel V{r), (59) 




where Ap X and A^ e denote the 4- vector potentials of the external and the self EM fields and z is a state to be 
suitably defined (see below). A clarification here is in order. The external EM 4-potential A^ xt ' (r) acting on the 
charged particle located at the 4-position r is assumed to be produced only by prescribed "external" sources, namely, 
excluding the particle itself, by the remaining possible EM sources belonging to the configuration space T r . Within 
the framework of special relativity, both the inertial term and the coupling term with the external field carry only 
local dependencies, in the sense that they depend explicitly only on the local 4-position r. They provide the classical 

(sel f) 

dynamics of charged particles in absence of any RR effect. On the other hand, the functional S 1 ^, associated to 
the EM self-interaction contains both local and non-local contributions. In particular, since the state z of a finite- 
size particle must include a 4-position vector r, it follows that S^ el ^ generally depends explicitly on two different 

4-positions, r and [r], to be properly defined (see below). The non-local property of Sq £1 ^ represents a characteristic 
feature of RR phenomena. 

From the relationship (|5"8"|) it follows that the Hamilton action functional for the treatment of the RR admits the 
decomposition (I12[) introduced by THM.l, namely it can be written as the sum of two terms, carrying respectively only 
local and both local and non-local dependencies. In order to prove that the same functional admits also a Lagrangian 
and a Hamiltonian representation in standard form it is sufficient to show that the self-coupling functional is symmetric 
in z and [z], in the sense defined in THM.l. For this purpose we need to determine explicitly the general expression 
of the 4-current and the self 4-potential for a rotating finite-size charged particle. 

The first step consists in constructing a covariant representation for the 4-current density. We follow the approach 
presented by Nodvik [20l ] . Thus, we consider an extended charged particle with charge and mass distributions having 
the same support dfl, to be identified with a smooth surface. Denoting by (s) the 4- vector position (with proper time 
s) of a reference point belonging to the internal open domain fl and by £ M a generic 4-vector of dQ, the displacement 
vector £^ is defined as: 



(61) 
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The particle model is prescribed by imposing the constraints of rigidity of dQ, namely for all <r M and (s) [201 ]: 



Z^Zii = const., 



(62) 
(63) 



where u^(s) = ^r M (s). In particular, we shall assume that mass and charge distributions are spherically symmetric 
and therefore characterized by a form factor / (£ 2 ) = / (£ M £ M ). This allows one to identify r M (s) as the center-point of 
dfl. The extended particle can in principle exhibit both translational and rotational degrees of freedom. In particular, 
the translational motion can be described in terms of (s). Instead, the rotational dynamics, which includes both 
space-time rotations associated to the so-called Thomas precession and pure spatial rotations, can be described in 
terms of the Euler angles a(s) = {(f(s), i?(s), i/>(s)}. It follows that, in this case, the Lagrangian state z must be 
identified with the set of variables z = (r a (s) ,a(s)). In view of these definitions it is immediate to prove that the 
4-current density for the finite-size particle can be written as follows: 



jW"(r)=?c 



where 



1 du a 
1 i—x 

as 



r" - r" (s) 



~u» v x v \ f{x 2 )5{x a u a ), 



(64) 



(65) 



and uj^ = lo^ u (s) is the antisymmetric angular velocity tensor [20], which depends on s through the Euler angles a(s). 
The term [l + d^"] contains the acceleration of r^ (s) and represents the contribution associated to the Thomas 
precession effect. This can be formally eliminated by using the properties of the Dirac-delta function, implying the 
identity: 



S(x a u a (s)) 



1 



ds 



6(s - si) 



1 



-S(s - s x ), 



ds 



where by definition s% = S\ (r) is the root of the algebraic equation 

u^O^-r" (sx)} = 0. 
As a result, the 4-current can be equivalently expressed as 



j(.«*f)»(r) =qc 



H-oo 



ds 



u^Sis - si) - -^ v x v 5(x a u a ) 



(66) 



(67) 



(68) 



The second step consists in constructing a Green-function representation for the EM self-potential A^ 1 ^^ in terms 
of the 4-current j( sel f)^(r). This technique is well-known. Thus, considering the Maxwell equations in flat space-time, 
in the Lorentz gauge A( sel f>@,g = 0, the self 4-potential must satisfy the wave equation 



|-|^(se;/)/j 



-in 



-J 



(self)ii 



(r), 



where □ represents the D'Alembertian operator and j( se/ /)^(r) is given by Eq. 



4-7T 

c 



d 4 r'G(r, r')j M/)Al (r'), 



(69) 

The formal solution of Eq. (|6"9"|) is 

(70) 



where G(r,r') is the retarded Green's function corresponding to the prescribed charge density. By construction, it 
follows that G(r,r') is symmetric with respect to r and r', and furthermore - since the particle is finite-size - both 
the 4-current and the self-potential are everywhere well-defined. 

From these general results, it is immediate to prove the following theorem. 

THM.3 - Symmetry properties of S^ eLf) (z, [z]) 

Given validity of Eq. i68\) for the covariant expression of the current density for a finite-size charged particle and of 
Eq. |70[) for the general expression of the corresponding EM self-potential, it follows that: 
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T3\) The functional Sq £1 ^ {z, [z]), defined in Eg. A60\) as an integral over the ^-volume element d 4 r, can be written 
as a line integral of the form 



/+oo 
dsL^{zM), (71) 
-oo 

where L^ el ^ represents the Lagrangian of the coupling with the EM self-field. This is defined as 



2 

4„ 



u M (5(s — si) oj^ x v 8(x a u a ) 

c 



f(x 2 ) A'G(r,/)jf . (72) 



T<?2,) The functional Sq £1 ^{z, [z]) contains both local and non-local dependencies in terms of the variational quan- 
tities z = z (s) and [z] = [z (s)] . In particular, it is symmetric in these local and non-local variables, in the sense 
stated in THM.l, namely 

S^ lf \z,[z]) = S^ lf \[z},z). (73) 

T3s) The functional Sq £1 ^ (z, [z]) contains at most only first-order derivatives of the variational functions z (s). 

Proof - T3i) The proof of the first statement follows by noting that the action integral Sq c1 ^\z, [z]) is a 4-scalar 

by definition. Hence, making explicit the expressions of A^ sel ^^ and jj? el ^ in Eq. (|60| according to the results in 
Eqs. fTPf and (fS5|) . by exchanging the order of the integrations and invoking the symmetry property of the Green 
function, the conclusion can be easily reached. In particular, the variational Lagrangian is found to be of the general 
form given in Eq. (|72p . 

T32) To prove the second statement we first notice that in Eq. ljTTI) both z and z' are integration variables, while 
by definition the variational quantities are identified with z (s) and [z (s)] = z 1 (s'). These dependencies are carried 
respectively by the charge current densities j( sei /)' I (r) and j( se '/)^(r'). The result is then reached by noting that the 
functional carrying the self-coupling terms is symmetric with respect to the integrated quantities, and in particular 
with respect to j^ sel f>^(^r) and j'( sei -W(r'). Hence, exchanging (z, j'( se ^)^(r)) with (V, j( se '/)^(r')) does not affect the 
form of the functional, with the consequence that Eq. (|T3")) is identically satisfied. 

T33) The proof of the statement is an immediate consequence of the representation for the current density j( sei /V(r) 
given in Eq. (p5|) . In fact, the term proportional to the acceleration in Eq. (|M|) and which is associated to the 
Thomas precession, does not appear in Eq. tpjS]) . thanks to the property of the Dirac-delta function indicated above 
in Eq.®. 

Q.E.D. 

An immediate consequence of THM.3 is that, thanks to THMs.l and 2, the variational treatment of the dynamics 
of finite-size charged particles subject to the EM RR effect admits both Lagrangian and Hamiltonian representations 
in standard form. In particular, in this case, it follows that the following identification must be introduced: 



L b =L^ elf \ (74) 



where is the Lagrangian defined above in THM.l. 



V. HAMILTONIAN THEORY FOR THE RR PROBLEM 

In this section, based on THMs.1-3 and the theory developed in Paper I, we proceed constructing the Hamiltonian 
formulation for the RR problem. For this purpose, it is convenient to recall the explicit form of the EM self 4-potential 
obtained in Paper I. For the sake of comparison with traditional approaches based on point particle models, here we 
also propose an alternative approach based on the Green function method. Remarkably, as pointed out in Appendix 
A, for the spherically-symmetric and non-rotating extended particle considered here, the self-potential is proved to be 
formally analogous to the well-known solution valid for point charges. This result holds, however, only in the external 
domain (with respect to dfl), where A^ el ^ (r) is found to admit the integral representation (see details in Appendix 
A): 

A^\r) = 2q f 2 dr'S(R a R a ). (75) 
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Here R a = r a — r a (s'), with r a and r' a = r a (s') denoting respectively the generic 4-position and the 4-position of 
the center of the charge distribution at proper time s' . As a fundamental consequence of the finite extension of the 
particle and the restrictions on the domain of validity of Eq . (f75j) . the resulting variational functional and Faraday 
tensor for the self-field turn out to be completely different from the point-particle treatment. In particular, the action 
integral becomes now a non-local functional with respect to the 4-position r. As pointed out in Paper I, this can be 
written as a line integral in terms of a variational Lagrangian L±(r, [r]) as follows: 

/+oo 
dfl£i(r,[r]). (76) 
-oo 

Here Ia(r, [r]) is defined as: 

Lr(r, [r]) = L M (r) + L% xt \r) + L^\r, [r]), (77) 

where 

L M (r,u) = -m c— — , (78) 
L { c Xt \r) = l^A^\r), (79) 

— (ext) 

are the local contributions respectively from the inertial and the external EM field coupling terms, with denoting 

the surface-averaged external EM potential (see Paper I for its definition). On the other hand, L^ el ^ represents the 
non-local contribution arising from the EM self-field coupling, which is provided by 

L { S elf \r, [r]) = f dr'MWR, - a% (80) 



c ds 



where the 4-scalar a 2 = is the radius of the surface distribution with respect to the center r M (s) and RJ 1 is 

defined as 

R a = r a (s)~r a (s'). (81) 

Notice that R a represents the displacement bi- vector between the actual position r a (s) of the charge center at proper 
time s and the retarded position r a (s') of the same point at the retarded proper time s' . It is immediate to verify 

that the representation of S^ el ^ in terms of L^ eL ^ given in Eq. (|80| satisfies the hypothesis of THM.l, and therefore 
the solution admits a Lagrangian representation in standard form. As already shown in Paper I and according to 
THM.l, this is obtained by setting 

L eff ee L M {r) + L^ xt \r) + 2L^ lf \r, [r]), (82) 

with Ljitf-(r), L^ xt \r) and L^ el ^ respectively given by Eqs. ([78]) - (j80 ]> . Then, the corresponding E-L equation is 
provided by the following covariant 4- vector, second-order delay- type ODE: 

du^s) q-j=(ext) ^dr u (s) q—{self) r , dr k (s) 

where 

«"00 = ^. (84) 
as 

Here the notation is as follows. Denoting by ee F^t^ + F/^f"^ the total Faraday tensor, Fjfu^ and F^ 1 ^ are 

respectively the "external" and "self Faraday tensors generated by A^ xt ^ and A^ el ^' , which carry the contributions 
due to the external sources with respect to the charged particle and the particle EM self-interaction. In particular, 

(ext) 

the 4-tensor F ' (r(s)) denotes the surface-average of the Faraday tensor associated to the external EM field, to be 
identified with 

F% xt) ^ d,Al ext) - d„A { ; xt \ (85) 
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(ext) 

with A v (r(s)) only generated by external sources with respect to the single-particle whose dynamics is described by 

(self) 

Eq. ([53"|) . Similarly, F^ k is the surface-average of the Faraday tensor contribution carried by the EM self 4-potential. 
In the parameter-free representation this is given by 



F ( ;: lf \r, W) = -Aqf \dr^JL5 ( R «R a - a 2 ) - dr' k ^S (iTR. a> 



(86) 



(self) 

As pointed out in Paper I, F^ k can also be parametrized in terms of the particle proper time s, by letting r = r (s) 

and [r] = r (s 1 ) in the previous equation, which also implies dr'^ = ds'-^r. This means that the non-locality in Eq. (|86[) 
can be interpreted as non-locality in the particle proper time. 

The remarkable feature of Eq. (|56T) is that the RR self- force (see the second term in the rhs of Eq. fSB")) ) contains non- 
local effects only through the retarded particle 4-position and not through the 4-velocity. This feature is fundamental 
for the subsequent fluid treatment, since it permits the evaluation in the standard way of the velocity moments, 
retaining the exact form of the RR self-interaction. 

The system of Eqs. ([83"l) and (|84|) defines a delay- type ODE problem of the form 

$ =X*(y,[r]) > 
y(s ) = yo, (87) 
y(s'o) = y s > , Vs e 7 S0 , S0 _ Sret , 

with so and s ret denoting respectively the initial particle proper time and the causal retarded proper time (see Paper 
I), and X/f the Hamiltonian vector field 

* H{ yM) = { dH '% PM \ - aH '>tf' [r]) }' <*> 

Denoting by y (s) = \ (yo, {y s ' i Vs € Is ,s —s ret } , s — so) the formal solution of the problem (|87j) . in the reminder 
we shall assume that the map 

yo^y(s) (89) 

is a diffeomeorphism of class C k , with k > 1. 

Based on these results, the Hamiltonian formulation is provided by the following theorem. 



THM.4 - Non-local variational and effective Hamiltonian functions for the non-rotating particle 

Given validity of THMs. 1-3, it follows that: 

T4i) The RR equation 1183]) for a non-rotating and spherically- symmetric charged particle admits the non-local 
Hamiltonian system {y = (r 11 ,p^),Hi}. Here p^ and Hi = Hi (r,p, [r]) are respectively the canonical momentum I2%jl 
defined with respect to the variational Lagrangian L x given in Eq. J77| ), and the corresponding non-local variational 
Hamiltonian 1121]) defined as the Legendre transformation of L\ . In particular, the variational non-local Hamiltonian 
121]) is identified with 

F lM ,W)-^^-^)(^-H, (90) 
where is the total EM ^-potential 

A,(r,[r]) = A { ; xt) (r)+3^ V) (r, [r]) , (91) 

(self) 

and from Eq. h80\) is the functional 

A ( ; elf) (r, [r]) = 2qj" dr'^R, - a 2 ). (92) 
T4-2) There exist P M and H e ff, defined respectively by Eas. \35]) and \34\ ), such that 

*«t k * M) - ^ ^ ~ l A ^) ( PM - H">) ' (93) 
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with A( e ff)n the non-local effective EM 4~ potential 



A 



{effhl (r,P)^A { ; xt) (r) + 2^r J >, [r]) 



(self) 



(94) 



(self) - 

ana A) L " defined in Eq.^M) 



T4.3) The effective and variational Hamiltonian junctions H e f f and Hi coincide when expressed in terms of the 
4-velocity d ; . 

Proof - The proof of T4i and T42 follows immediately by applying THMs.l and 2 with the variational Lagrangian 
L\ given by Eq. ([77|) . In particular, this yields 



ds 



(ext) 



and 



m c- 



ds 



A 



—r(ext) 



.4 



(95) 



(96) 



The corresponding Legendre transformations then provide respectively Eq. ([9H)l and Eq. (j9"3")) . Finally, by direct sub- 
stitution of Eq. (j9"5"l) into Eq.(j90]) and Eq.([96]) into Eq. (l93"l) . one obtains that 



H, 



eff 



#1 



m c dr^(s) dr fJ '(s) 
2 ds ds 



(97) 



which proves also the last statement. 
Q.E.D. 



We remark that the Hamilton equation in standard form expressed in terms of H e f t and P M are differential equations 
of delay- type, as a consequence of the non-local dependencies appearing in H e f t which are characteristic of the RR 
phenomenon. In this case, for the well-posedness of the solution the initial conditions in the interval I = [so — s r et, so] 
must be defined, with so the initial proper time and s re t a suitable retarded time. However, if the assumption of 
inertial motion in the proper time interval Iq — [—00, sq] holds, then the mapping 



yo = y (so) -> y (s) = T SOtS y , 



with y = (r M , Pu), defines a classical dynamical system (see Paper I), and this dynamical system is Hamiltonian. 



VI. A HAMILTONIAN ASYMPTOTIC APPROXIMATION FOR THE RR EQUATION 

In this section we intend to carry out in detail a comparison of the present approach for extended particles with 
the customary point-particle treatments leading to the LAD and LL equations. For this purpose, asymptotic approx- 
imations of the exact RR self-force ([55]) are investigated. 

The issue has been partially discussed in Paper I. As pointed out therein, an asymptotic approximation of the exact 
RR equation (|83[) can be obtained in validity of the short delay-time ordering, namely requiring 



< e 



Sret 



< 1, (99) 



where s re t = s — s' , with s and s' denoting respectively the present and retarded particle proper times. This permits 
two different possible strategies, respectively based on Taylor expansions performed with respect to s (present-time 
expansion) or s' (retarded-time expansion). In particular, adopting the present-time expansion for the RR self-force 
(|86|) . the delay- type ODE ([831 can be reduced, in principle, to an infinite-order differential equation. Instead, by 
truncating the Taylor expansion to first-order in e, ignoring mass-renormalization terms and taking the point-particle 
limit a — > 0, in this way the customary expression for the LAD equation is recovered (see THM.3 of Paper I). 

As remarked in the Introduction, the resulting asymptotic approximation (given by the LAD equation) is non- 
variational and therefore non-Hamiltonian. In addition, contrary to the exact RR equation obtained here, the LAD 
equation, as well as the related LL approximation, both fail in the transient time intervals occurring when the external 
EM field acting on the particle is turned on and off. To elucidate this point, let us consider the dynamics of a charged 
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particle which is in inertial motion in the past for all s < sq and from s = sq is subject to the action of an external 
EM field. Then, by construction, it is immediate to show that in the transient time interval Iq — [sq, Sq + s ret ] the 
exact RR self-force (|83|) is manifestly identically zero. In fact, in the case of inertial motion in the past (namely 
U(i(s') — const.) the RR self-force vanishes in such a time interval (see THM.l in Paper I). In contrast, both the LAD 
and LL equations predict incorrectly a non- vanishing RR self- force. The same kind of inconsistency (for the LAD and 
LL equations) arises when the analogous transient time interval corresponding to the turning-off of the external EM 
field is considered Q- 

Therefore, the issue arises whether an alternative asymptotic approximation can be determined (for the exact RR 
equation) which simultaneously: 

1) overcomes this deficiency, by taking into account consistently relativistic finite delay-time effects characteristic 
of the RR phenomenon; 

2) is variational and admits a standard Hamiltonian formulation. 

In this section we propose a solution to this problem, by performing a retarded-time expansion, which provides an 
alternative to the LAD and LL equations. 



A. The Hamiltonian approximation 



For definiteness, let us assume that the external force acting on the particle is non- vanishing only in a finite proper- 
time interval / = [so,Si]. Then, in validity of the ordering (|99[) . we require that the external EM force is slowly 
varying in the sense that, denoting r' = r' 1 (s') and r = r^ (s), 

^TV)-^fV) ~0(e), (100) 
(F ( »f(r')-F ( ;f{r)) h -0(e), (101) 

(F^ (r') -F^ (r)) ' ~ 0(e). (102) 

V / ,hk 

Then, the retarded-time Hamiltonian approximation of the RR equation is obtained by performing a Taylor expan- 
sion in a neighborhood of s' . The result is summarized by the following theorem. 

THM.5 - First-order, short delay-time Hamiltonian approximation (retarded-time expansion). 

Given validity of the asymptotic ordering $9Sfy and the smoothness assumptions MOty - HU^I for the external EM 
field, neglecting corrections of order e™, with n > 1 (first- order approximation), the following results hold: 
T5\) The vector field 

q —(self) , dr k (s) 

G n=~ F fl k (r(s),r{s)) ^ (103) 

appearing in Eq. &83\) can be approximated in a neighborhood of s' as 

9n (r (a')) = {-^oEmc^u^ (s') + ^ (r (s'))\ , (104) 
to be referred to as retarded-time Hamiltonian approximation, in which the first term on the rhs identifies a retarded 

2 

mass- correction term, m EM = denoting the leading-order EM mass. Finally, g'^ is the J^-vector 



4. (>•(»'» =4t 



d 2 d 2 

(s') - u ^s')u k {s')—u k {s 1 ) 



(105) 



T52) The corresponding RR equation, obtained replacing G M with the asymptotic approximation g^ is vari- 

ational, Lagrangian and admits a standard Lagrangian form. Let us denote with r' a = rg (s') the extremal particle 
world-line at the retarded proper time s' . Then, in this approximation the corresponding asymptotic variational La- 
grangian and effective Lagrangian functions coincide. Both are defined in terms of the asymptotic approximation 
L^ e ^y m (r,r' ), replacing L^ cl ^ . To leading-order in e, this is found to be 

L [ c: l alUry n )=g,{r',)r». (106) 
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T5s ) The asymptotic approximation given by Eq. {10$ is also Hamiltonian. The asymptotic variational and effective 
Hamiltonian functions coincide and are given by 

dr^ 

with 

L hasym (ry ) = L M (r)+L% xt \r) + L ( ^ ym (r,r> ), (108) 

and now 

dL hasym 

u ds 

Proof - T5i) The proof can be carried out starting from Eq. ([83|) and performing explicitly the Taylor expansion in 
a neighborhood of s' = s — s re t- For a generic analytic function / (s), this yields the power series of the form 

k=0 

In particular, for the 4- vectors dr ^ s ^ and R k one obtains respectively the asymptotic approximations 

dr,(s) „ dr, (s') d 2 r, (s') (s s') 2 d 3 r, (s') 

ds ~ ds' +[ ' ds' 2 + 2 ds' 3 + U ^> 

and 

~. , .^dr k (s') (s-s') 2 d , „ (s — s') 3 d 2 . , ,. „ , As , 
R k = (s — s')-^ + L_^_ u * (*') + { -^^u* (s') + O {^) , (112) 

while for the time delay s — s' = s re t the leading-order expression 

s-s' = a + 0(e 2 ) (113) 



holds. By substituting these expansions in Eq. ([55|) . the asymptotic solution given by Eq. (ll04[) can be recovered. 

T52)-T53) The proof follows by first noting that T^, el J^ ym contributes to the Euler-Lagrange equations only in terms 
of the local dependence in terms of r. Then, in this approximation the canonical momentum becomes 

dr u (s) q—iext) , s _ /, , «\ 

while the asymptotic Hamiltonian reduces to 

H hasym (r,p, r' ) = -L- (p„ - ^A { ; xt \r)) (> - U { ^(r)) + g, (r' ) r". (115) 

The corresponding Lagrangian and Hamiltonian equations manifestly coincide with Eq. (|83l) once the approximation 
f)104|) is invoked for the vector field G M . 
Q.E.D. 



B. Discussion and comparisons with point-particle treatments 



The asymptotic Hamiltonian approximation, here pointed out for the first time (see THM.5), preserves the basic 
physical properties of the exact RR force (|83|) . In fact in both cases, the RR force: 

1) is non-local, depending on the past history of the finite-size charged particle; 

2) admits a variational formulation; 

3) is both Lagrangian and Hamiltonian; 
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4) satisfies the Einstein Causality Principle and, when applicable, the Newton Principle of Determinacy (see also 
Paper I); 

5) describes correctly the transient time intervals in which the external force is turned on and off (sudden force). 

For these reasons, physical comparisons based on the retarded-time Hamiltonian asymptotic approximation are 
meaningful. In particular, here we remark that the present approach departs in several ways with respect to point- 
particle treatments based on the LAD and LL equations. More precisely: 

1) The same type of asymptotic ordering is imposed, which is based on the short delay-time ordering (|99[) . However, 
in contrast with the LAD and LL equations, the expansion adopted in THM.5 and leading to the retarded-time 
Hamiltonian approximation can only be performed based on the knowledge of the exact RR force for finite-size 
particles. 

2) Unlike the LAD and LL equations, the asymptotic Hamiltonian approximation carries the information of the past 
dynamical history of the charged particle through the retarded time s' . Therefore, the dynamical equation written 
adopting the approximation (I104[) is still a delay-type second-order ODE. The construction of its general solution 
becomes trivial in this case, since the self- force is considered as an explicit source term evaluated at proper time s' . 

3) The asymptotic approximation provided by Eq. (|104p cannot be regarded as a point-particle limit. In fact, the 
retarded mass-correction term would diverge in this limit. 

4) The exact RR equation satisfies identically by construction the kinematic constraint u^u^ = 1. The same 
constraint is satisfied to leading-order in e also both by the retarded and present-time asymptotic expansions (and 
hence also the LAD equation). 

5) The variational principle introduced in THM.5 is subject to the constraint that the past history is considered 
prescribed in terms of the extremal world-line. This requirement is consistent with the initial conditions for the 
RR equation, which is a delay-type ODE depending only on the past history of the particle. This requires that the 
world-line trajectory is prescribed in the past, namely in the time interval / = [— oo, So]. Since, however, the initial 
proper time So is arbitrary, it follows that r (s) can be considered prescribed also in the time interval I' = [— oo, s 1 ]. 
In particular, if for all s < so the motion is assumed to be inertial, the initial-value problem associated to the RR 
equation written in terms of the retarded asymptotic self-force (|104p is well-posed, in the sense of the standard 
Newton Principle of Determinism, as discussed in Paper I (see in particular THM.4 presented there and dealing with 
the existence and uniqueness of solutions for the exact RR equation). 

6) One might think that the same type of constrained variational principle, of the kind adopted in THM.5, could 
be adopted also for the exact RR equation. However, this belief is wrong. In fact, since the variational functional ([77)1 
is symmetric with respect to the local and non-local world-line trajectories, there is no distinction between past and 
future. Since future cannot be prescribed, such a constrained variational principle for the exact equation is forbidden. 
On the contrary, the extremal RR equation (|83p is obtained by imposing also the Einstein Causality Principle, and 
therefore it depends only on the past history. 

7) Despite some formal similarities between the retarded-time Hamiltonian approximation versus the corresponding 
LAD and LL equations, the latter cannot be recovered even in the framework of some kind of constrained variational 
principle. In fact this would require to consider prescribed for example, second or higher-order proper-time derivatives 
of the particle position vector (namely the acceleration and its derivatives). This viewpoint is manifestly unacceptable, 
because it would amount to constraint the present state of the particle at proper time s. 

8) The previous argument justifies, in turn, the introduction of the short delay-time asymptotic approximation 
given in THM.5. This is performed directly on the RR force, namely the 4-vector G M entering the RR equation itself. 
In this way the variational character of the RR problem is preserved. It follows that the corresponding variational 
functional as well as the Lagrangian and Hamiltonian functions for the asymptotic RR equation are constructed only 
"a posteriori" . 

9) Another advantage of the new representation (I104p with respect to the customary LAD and LL equations is that 
it permits the approximate treatment of the solution also in the transient time intervals after the turning-on or the 
turning-off of the external EM field. In particular, in contrast to the LAD and LL equations, it predicts a vanishing 
RR self-force in the turning-on transient phase Iq — [so, sq + s ret \. 

10) Finally, it should be remarked that the retarded asymptotic self-force (|104p cannot be trivially obtained from 
the corresponding local asymptotic representation performed at proper time s and leading to the LAD equation by 
simply exchanging s with s' (or by a further Taylor expansion). Indeed, the relationship between the two can only 
be established based on the exact form of the self-force. 
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VII. COLLISIONLESS RELATIVISTIC KINETIC THEORY FOR THE EM RR EFFECT - CANONICAL 

FORMALISM 

In this section we proceed with the construction of the relativistic classical statistical mechanics (CSM) for a 
collisionless plasma with the inclusion of the EM RR effect. In particular we shall prove that the mathematical 
formalism introduced in the previous sections to deal with symmetric non-local interactions allows one to obtain a 
convenient formulation for the kinetic theory describing such a system and for the corresponding fluid representation. 
The derivation is based on the property of a symmetric non-local system represented by a hnite-size charged particle 
of being Hamiltonian with respect to and H e f f . 

In view of the peculiar features of the non-local RR phenomenon and the related delay-type differential Hamiltonian 
equations, it is instructive here to adopt an axiomatic formulation of the CSM for relativistic systems with the inclusion 
of such an effect. We shall assume that the latter are represented by a system of classical finite-size charged particles 
subject only to the action of a mean- field external EM force and a non-local self-interaction. We intend to show that, 
using the Hamiltonian representation in standard form given above, the explicit form of the relativistic Vlasov kinetic 
equation can be obtained for the kinetic distribution function describing the statistical dynamics of such a system. 
Therefore, the problem is reduced to a Vlasov-Maxwell description for a continuous distribution of relativistic charged 
particles. 

For dcfinitcness, let us consider the non-local Hamiltonian dynamical system in standard form {y, H e f /} given above. 
This is characterized by the superabundant state vector y = (r^,P^) spanning the extended 8th-dimensional phase- 
space r and with essential state variables yi (y) spanning the 6th-dimensional reduced phase-space T\. Introducing 
the global proper time s, T\ (s) is defined as 



where |w| = y / u a u a and s (y) is the world- line proper time uniquely associated to any y. By assumption, Ti (s) is an 
invariant set, i.e., Ti (s) = Ti for any s e M. Next, let us consider the Hamiltonian flow T S(uS defined in Eq. |98|) . By 
construction the dynamical system is autonomous, namely the flow is of the form 



The existence of the dynamical system T S0)S for the state y (s) has been proved in Paper I. This requires that in 
the proper time interval Iq = [—00, sq] the motion of each charged particle is inertial, namely the external EM field 
vanishes in the same interval. As a result of Eq. (|98[) . any point in the phase-space T spanned by y or yo is associated 
to a unique phase-space trajectory, namely such that y = y(s), for any y G T. Due to (|98|) there exists necessarily 
yo = y (so) which is mapped in y (s) . Viceversa, for any s€l there exists a unique y = y (s). However, we notice 
here that for the axiomatic formulation of the CSM for the RR problem the assumption of existence of the dynamical 
system T So , s is not a necessary condition. In fact, it is immediate to prove that the minimal requirement is actually 
provided only by the existence of the diffeomeorphism (|89p defined above. 

Now, for a prescribed So £ K let us consider the set B (sq) ^ Ti, with B (So) an ensemble of states yo, each one 
prescribed at the initial proper time sq = s'a- Its image generated at any s = s e I by the flow T S0[S , for each 
trajectory, is 



where s and sq denote now the global proper times Sand s~o> 
We introduce the following axioms. 
AXIOM #1: Probability on K (T). 

Let K (Ti) be a family of subsets of Ti which are L-measurable. We define the probability of B (s) e K (Ti) as the 



r i («) = {y : y e r, \u\ = 1, s (y) = s, ds (y) = y/g^ydr^dr^ , 



(116) 



T So , s yo = x(yo,s - s ) ■ 



(117) 



B = B (s) = T S0 . S B (s ) , 



(118) 



function 



P(B) -.KiT^^ [0,1 



(119) 



such that it satisfies the constraints 



P(0) 



1, 
0. 



(120) 
(121) 



OO 




(122) 
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with {B. L £ K (Ti) , i £ N} being an arbitrary family of separate sets of K (Yi). 
AXIOM #2: Probability density. 

For any B (s) £ K (T{) and for any state y = (r^, P^) there exists a unique probability density p (y) > on Ti such 
that 

P (B (s)) = J dyp (y) S(\u\ -l)S(s-s (y)) S B(s) (y) , (123) 

where dy — dr^dP^ is the canonical measure on T and Sb( s ) (y) is the characteristic function of B (s). Furthermore, 
s (y) is a particle world-line proper time, while s = So + As, with As an invariant proper time interval independent 
of So- We notice that s (y) can be equivalently parametrized in terms of the observer's coordinate time r , namely: 

ds(y)^dr^g^ — — . (124) 

AXIOM #3: Equiprobability. 

Then, the equiprobability condition requires that, for all B (sq) and for all s, Sq £ I C K, 

P(B(s))=P(B(s )). (125) 



We remark that in the integral (|123[) the two Dirac-delta functions can be interpreted as physical realizability 
conditions, required to reduce the dimension of the volume element dy defined on the extended phase-space V. 
We can now introduce the following theorem, concerning the validity of the Liouville equation for p(y). 

THM.6 - Relativistic Liouville equation for p(y). 

Given a Hamiltonian system {y, H e ff} and imposing the validity of Axioms #l-#3, it follows that the probability 
density p(y(s)) is a constant of motion, namely for any s,so £ M (to be intended now as world-line proper times) 
and for any yo £ T 



p(y(s)) = p(yo) , 

to be referred to as the integral Liouville equation. This can also be written equivalently as 

d 



ds 



p(y(*)) = o, 



(126) 



(127) 



to be referred to as the differential Liouville equation. As a consequence, introducing the kinetic distribution function 
(KDF)f(y) 



f(y) = p(y)N, 



(128) 



with N being the total number of particles in the configuration space of B C K(T), it follows that also f (y) satisfies 
the Liouville equation 

Proof - We first notice that, from Axiom #1, by changing the integration variables we can write Eq. (ll23p as 



(129) 



P (B (s)) = / dyp (y) S(\u\ - 1) 5 (s - s (y)) S B(s) (y) = 

dy (s) 



dyo 



dyo 



p (y (s)) 5 (\u\ - 1) 5 (s - s (y)) 5 B(so) (y (s )) 



with 



dyo 



being the Jacobian of the variable transformation from y (s) to yo- On the other hand, since the system 



{y,H e ff} is Hamiltonian, it follows identically that 



9y(s) 



<9yo 



= 1. Hence, invoking Axiom #2 we can write 



L 



dyo [P (y (s)) 5 (\u\ -l)5(a-8 (y)) - p (y„) S (Kl -l)S(s -s (yo))] 6 B{so) (y (s )) = 0, 



from which it must be that 

p (y (s)) S (\u\ ~l)S(s-s (y)) = p (y ) 6 (\u \ -l)S(s -s (y )) . 



(130) 



(131) 
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On the other hand, by construction it follows that 



*(M-i) 



1 


d\u 




d\u \ 


1 


-6 


ds 
ds 





*(|uo|-l) = *(|uo|-l) ; 



S{s - s (y )) =6(s -s (j/ )) 



In fact, by definition the 4- velocity is normalized to 1 at all proper times, so that 



d\u\ 



d\u \ 



with As being independent of the initial value sq, and hence 



1 too. 



Finally, because of these conclusions, from Eq. (|13ip it follows that 

p(y(s)) = p(yo) , 



(132) 
(133) 

= 1. Furthermore, s = sq + As, 

(134) 



which represents the Liouville equation in integral form. By differentiating with respect to s the equivalent differential 
representation follows at once. An analogous equation holds manifestly also for the KDF / (y). 
Q.E.D. 

We conclude noting that, formally, the Liouville equation for non-local Hamiltonian systems in standard form is 
analogous to that characterizing local Hamiltonian systems. Such an equation can be viewed as a Vlasov equation for 
a relativistic collisionless plasma, in which each particle is subject only to the action of a mean-field EM interaction, 
generated respectively by the external and the self EM Faraday tensors. By definition, in this treatment the latter do 
not include retarded binary EM interactions. It follows that, in terms of the Lagrangian equation (|127p . the probability 
density p (y (s)) is parametrized in terms of the single-particle phase-space trajectory {y (s) , s £ I}. Hence, it advances 
in (proper) time s by means of the canonical state y (s) as determined by the Hamiltonian equations of motion (|87l) . 



A. Vlasov-Maxwell description 

To define a well-posed problem, the relativistic Vlasov equation (|127[) must be coupled to the Maxwell equations, 
which determine the total EM field produced by all the relevant sources. In particular, in order to determine the 
external Faraday tensor F^t xt \ the corresponding EM 4-potential A^ 3 *' must be determined. In the Lorentz gauge, 
this is prescribed requiring it to be a solution of the Maxwell equations 

Arrr 

where j( ext )^(r) is identified with the total current density 

j( ext )»(r) ee qj d A u5 (M - 1) vTf (y) + j( cMs )^( r ). (136) 

Here, the first term is the Vlasov 4-current density, namely the velocity moment of / (y) carrying the non-local phase- 
space contributions which yield the collective field produced by the plasma. The second term, instead, is produced 
by possible prescribed sources located outside the plasma domain. Therefore, in the Vlasov-Maxwell description the 
total EM 4-potential acting on a single particle must be considered as represented by A v — A^ xt ^ + Ai? el ^ , where 
Ai self) is given by Eq.(75) and A ( u ext) is the solution of Eq.([l35]). 

Therefore, the dynamical evolution of the KDF along a single-particle phase-space trajectory depends both explicitly, 
via A^ el ^\ and implicitly, via the 4-current j( ea;t )^(r), on the whole Faraday tensor ee Fj^ 1 ' + F^t el ■ In this 
way contributions which are non-local both in configuration and phase-space are consistently included in the theory. 
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VIII. FLUID MOMENT EQUATIONS 



We now proceed to compute explicitly the relativistic fluid moment equations which follow from the Liouville 
equation. To this aim, the relativistic Liouville equation is conveniently written as a PDE (Eulerian form) 



where G M (y) is defined by Eg. (1551) . or as an ODE (Lagrangian form): 

dr»df(y(s)) , du^df(y(s)) _ 



ds dr^ 



ds du b 



0, 



with y (s) being the phase-space trajectory of a particle. Then, the relativistic fluid equations related to the 
equation are defined as the following integrals over the momentum space: 



d 4 uS(\u\ -1)G 



U »d_m + G » {y) df(y) 



du„ 



0. 



Similarly, the corresponding fluid fields are defined as 

d l u5(\u\-l)Gf (y), 

with G = 1, u^, u^u v , ... and u M is the 4- velocity. In particular, we shall denote 

n(r) = J d 4 M ,5(M-l)/(y), 
N»(r) = n (r) W* (r) = J d A u8 {\u\ - 1) u M f (y) , 
T^(r) = / d 4 uS(\u\-l)u^f(y), 



(137) 

(138) 
Liouville 

(139) 
(140) 

(141) 
(142) 
(143) 



to be referred to as the number density, the 4~flow and the stress- energy tensor. 
It is immediate to prove that the corresponding moment equations are as follows. 



Continuity equation 

For G = 1 the Liouville equation provides the continuity equation 

fyJVT" (r) = 0. 

Energy-momentum equation 

For G = u u the Liouville equation provides the energy-momentum equation 

d ll T^{r)=F^ ot) (r)N,(r), 



(144) 



(145) 



where, from Eq. (f83j) we have that 



F v » ( r ) = 1 
( to *) 1 ' ~ m c 2 



(146) 



is the total EM force, with F^ sel ^ u ^ containing the retarded non-local contributions arising from the EM RR effect. 



We remark the following properties. 

1) As a consequence of the Hamiltonian formulation in standard form, the fluid equations obtained from the kinetic 
equation with the inclusion of the RR effect are formally the same as in the usual treatment for local systems. 

2) The contribution of the RR effect to the fluid equations is contained explicitly in the source term in the rhs of 
Eq. (|145[) , and also implicitly in the definition of the fluid fields. In fact, by assumption, the KDF is a function of 
the effective Hamiltonian state y = (r M ,P M ), which depends on the retarded self-potential. Hence, the fluid fields 
defined by Kqs. (|141| )- (|143[) must be interpreted as the fluid fields of the plasma which is emitting self-radiation and is 
therefore subject to the RR effect. 
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A. The implicit contribution of the RR self-force 



It is worth discussing the features of the theory in connection with the implicit contribution of the RR effect 
contained in the definition of the fluid fields. In particular, here we show that such contribution can be made explicit 
and an analytical asymptotic estimation of it can be give provide some suitable assumptions are imposed on the 
physical system. This concerns the case in which the contribution of the self-potential is small in comparison with 
the external EM potential in the KDF. In these circumstances, the exact KDF can be Taylor expanded as follows: 



/ (y) ~ / (y„ c ) + (y 



df(y) 



(147) 



y=y„c 



where y nc = (r^,^) is the state which is canonical in absence of the EM self-field. It is clear that, by construction, 
only the canonical momenta are involved in this expansion, since the configuration state is left unchanged by the 
presence of the self-force. Therefore, from the form of the previous expansion it follows that the first term of the 
series, namely / (y n c), does not contain any contribution from the RR self-field. Consider, for simplicity, the Taylor 
series to first order. Then, the corresponding fluid fields can be decomposed as follows: 



n (r) 
N» (r) 



no (r) + 
(r) 



ni (r) 



1 Q 



' 1 1 



(r), 
(r), 



(148) 
(149) 
(150) 



where 



n (r) 
ni (r) 



d 4 U <5(M-l)/(y„ c ), 
d\5 (\u\- I) (y-y nc ) 



df(y) 



Qy 



y=y« 



d 4 uS{\u\ - 1) 



df(y) 



=IH> 



(151) 
(152) 



and similar definitions hold for the other two fluid fields. 

To illustrate the procedure, let us consider, for example, the case of a relativistic Maxwellian distribution of the 
form ED 



Sm (y) = 



l 



(27170' 



■ exp 



T 



(153) 



where /i, P M , and T are respectively the chemical potential, the canonical momentum and the fluid 4- velocity and 
temperature. Then, in terms of the previous expansion, we obtain for the density 



n-o (r) 



nx (r) 



4irm 2 cT 



{2ttK) 



3-K2 [ -ppr I ex P 



T 



—r(ext) 
\ 

T c T 



{self) 

«o (r) 



2q A)7~ J 'U^ 



T 



(154) 



(155) 



with K2 (^f-^J being the modified Bessel function of the second kind. As can be seen, the effect of the RR self-field 



appears only in n\ (r) through the integral over the non-local dependencies contained in the potential 
follows that for a Maxwellian KDF the 4-flow (r) can be written as 

N» (r) ~ [n (r) + n x (r)] (r) , 

while the expansion terms of the stress-energy tensor T^ v (r) are given by 

T^{r) = \n eU^U" - poA^ , 
T? v (r) = \n x eU»U v - px^ v . 



(self) 



It 



(156) 

(157) 
(158) 
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Here the notation is as in Ref . [2lj . Thus, A A " y is the projector operator A Ml/ = rf v — c 2 U tl U v , e is the energy per 
particle 




(159) 



and from the definition of the pressure as p — nT it follows that 

Po(r) = n (r)T, (160) 

Pl (r) = n 1 (r)T = -^A i ° elf) U»no(r). (161) 
c p 

Finally, let us consider how the fluid equations are modified from the introduction of the series expansion (|147l) . 
Substituting the relations (jl48|) - (|150[) into the moment equations, for the continuity equation we get 

(r) = -d^N? (r) , (162) 

and for the momentum equation 

wr^p&^-wr- (163) 

In this way, on the lhs we have isolated the terms of the "unperturbed fluid" , namely the physical observables 
corresponding to a charged fluid in absence of RR. On the other hand, the asymptotic contributions of the RR effect 
have been isolated on the rhs, which allows one to interpret them as source terms due to extra forces acting on 
the unperturbed fluid. In particular, the presence of the RR acts like a non-conservative collisional operator, if we 
interpret it as a sort of retarded scattering of the fluid (and therefore, of the single particles at the kinetic level) with 
itself. 



IX. LAGRANGIAN FORMULATION OF THE FLUID EQUATIONS 

An important issue concerns the treatment of the non-local contributions appearing in the fluid equations both in 
the definitions of the fluid fields and in the source term in the momentum equation. This requires, in particular, the 

(self) (self) 

explicit representation of the self-potential and the EM self- force defined respectively in Eqs.(|92p and 



(|86|) . In fact, in the previous sections these non-local contributions have been written in a parameter- free representation 
(integral form), so that they do not depend on the retarded particle velocity. This allowed us to perform the velocity 
integrals in a straightforward way, only in terms of local 4- velocities, in agreement with the formalism adopted for 
the Hamiltonian formulation in standard form. 

To treat these non-local terms it is first convenient to represent the fluid moment equations in Lagrangian form, 
describing the dynamics of fluid elements along their respective Lagrangian path (LP). By substituting the definition 
(|142| in Eq. (|144p we obtain the corresponding Lagrangian form of the continuity equation, given by 

—n + nd^ = 0, (164) 

where -jy- = (r (s)) is the convective Lagrangian derivative along the LP of the fluid element parametrized in 

terms of the arc-length s, and U 1 * (r (s)) = dr d ^ ■ Similarly, writing the stress-energy tensor T^ v (r) as T^ v (r) = 
nU^U 1 ' + P^ v (r), with P A " y (r) = T^ v (r) — nU^U v ', the energy-momentum equation (|145l) can be represented in 
Lagrangian form as follows: 



"7^'" "l-"Zt i: ill- ci or,) 



Analogous results can be given for the asymptotic equations (|162[) and (|163|) . 

With the introduction of the LPs, the paramctrization of the non-local contributions can be easily reached in terms 

(self) 

of the LP arc-length s. Consider, for example, the self-potential . This can be expressed as 

A ( ; elf) (r, [r]) ee 2qj' ds' d -^8{R»R, a 2 ), (166) 
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dr' 



where by definition now = (r (s')) is defined along a fluid element LP. Then, by expressing the Dirac-delta 
function as 



*(JJ"iZp - a 2 ) 



2i? a f/ a 



5 (s' - S + S ret ) ■ 



(167) 



r(self) 



it follows that A can be equivalently written in the integrated form as 



R a U a (r(s')) 



(168) 



s'=s-s ret 



with R a being the displacement vector defined along a LP. In particular, in agreement with the Einstein Causality 
Principle, the retarded time s re t = s — s' is the positive root of the delay-time equation 



R^Rf, - cr 2 = 0. 

An analogous derivation can be carried out also for the self-force F „ fc , giving the following result 



(169) 



F { X S) (r, [r]) = -2«- 



1 D 



R a U a {s') 



Ds 



(r(s')) 



(170) 



s'=s— s re t 



where 



U^(r{s'))R k -U k (r [s'))Ru 



R a U a {r(s')) 

Again, this expression must be intended as a parametrization defined along a fluid element LP. 



(171) 



We conclude by commenting on the following remarkable aspects of the theory presented here. 

1) The fluid equations with the inclusion of the non-local effect related to the EM RR have been derived in a closed 
analytical form in both Eulerian and Lagrangian formulations. In particular, it follows that the fluid dynamics of the 
non-local kinetic system is intrinsically non-local too. 

2) Non-local contributions of the RR appear both in explicit and implicit contributions, through the definitions of 
the fluid fields as velocity moments of the KDF. 

3) From the point of view of the fluid description, it follows that the natural setting for the treatment of the non- 
local fluid equations is given by the Lagrangian formulation and the concept of LPs. This is a consequence of the fact 
that the exact moment equations are of delay-type. In fact, in order to properly deal with the non-local contributions 
of the RR the parametrization of the retarded effects in terms of the arc-length of the corresponding LPs is needed. 
It follows that the dynamics of a generic fluid element along its LP is related to the EM RR effect produced at the 
retarded time along the LP itself. 



X. ASYMPTOTIC APPROXIMATION 



In the previous sections we derived an exact formulation for both kinetic and fluid theories describing systems 
of relativistic charged particles subject to the EM RR self-interaction. In particular, we have pointed out that the 
kinetic and fluid equations are of delay-type, and therefore intrinsically non-local, due to the characteristic feature of 
the RR effect of being a non-local retarded effect. The retarded proper time is determined by Eq. (|169j) in agreement 
with the causality principle. Notice that this equation has formally the same expression for the single-particle or the 
kinetic dynamics and for the fluid equations in Lagrangian form (see also Paper I) . By inspecting Eq. (|169[) it is easy 
to realize that the order of magnitude of the delay-time is approximately s re t ~ c/c, and therefore very small for 
classical elementary particles. The smallness of the retarded time may represent a serious problem for the practical 
implementation of the exact theory presented here. In fact, the retarded time associated to the RR can be orders 
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of magnitude smaller than any other characteristic time for most of relevant physical situations. The question is of 
primary importance, for example, for the actual numerical integration of the exact fluid equations. 

In view of these considerations, in this section we provide asymptotic estimations of the non-local terms appearing 
in the moment equations, which allow one to overcome the difficulty connected with the finite delay-time intervals 
carried by the RR phenomenon. This requires to introduce a suitable asymptotic expansion of the exact non-local 
terms by means of approximations in which the self-interaction contributions are all expressed only through local 
quantities. The result has potential interest also in relation to the use of Eulerian integration schemes for the fluid 
equations with the inclusion of the RR effect. 

Specifically, the present analysis requires to develop an asymptotic approximation which involves the treatment 
of the delay-time s ret . This is accomplished within the short delay-time ordering approximation given by Eq.(|99p. 
In the following we shall work adopting the Lagrangian representation form for the fluid equations. To perform the 
asymptotic expansion, we assume that both the external EM field acting on each fluid element and the macroscopic 
fluid fields associated to the kinetic system are smooth function of the coordinate 4-position vector r a , namely they are 
of class C k , with k > 2. The result of the asymptotic approximation for the terms associated to the RR self-interaction 
is provided by the following theorem. 

THM.7 - First-order, short delay-time asymptotic approximation (present-time expansion). 

Given validity of the asymptotic ordering A99\) and the smoothness assumptions for the external EM and the fluid 
fields, neglecting corrections of order e", with n> 1 (first-order approximation), it follows that: 

T7i ) The retarded self-potential A ^ defined in Eq. U68\) can be expanded in a neighborhood of s as follows: 



-j(self) _ -j(self) 
A fi — A n 



where the present-time leading-order contribution 

(self) 



(self) 



1 + 0(e)] 
is given by 



A. 



\u,(r{s))-^-U,(r{s)) 



(172) 



(173) 



with being the convective derivative along a fluid element Lagrangian path. 
TI2) Concerning Eg. A165\ ), let us define the vector field as follows: 



K„ 



m c 



q_—( se if) 

2 r iiv u ! 



(174) 



with ^ defined in Eq. jl70\ ). Then, in a neighborhood of s, can be expanded as follows: 

K^= *g fl [l + 0(e)], 
where the present-time leading-order contribution is given by 

1 q 2 D 



1 a m c^ Ds 



U^{r{s))+g^ 



(175) 



(176) 



with g^ denoting the ^-vector 



9ix = o 



2 , 2 



3 m D c 2 



~U,-U^ S )U k { S )^U k 



Ds 



Ds 2 



(177) 



Proof - The proof of T7i) and T72) can be reached by introducing a Taylor expansion in terms of the retarded 
time s' for the relevant quantities appearing in Eqs. (|168p and (1170[) . In particular, for the 4- velocity U ^ (r (s')) and 
the displacement vector R k we obtain respectively 



D 



(s-s') 2 D 



U„ (r (s>)) - !/„ (r (a)) -{a- (r («)) + ^^j^I U » ( r + (^) 



and 
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:U k + e 4 , 



(178) 
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while for the time delay s — s' = s ret we get 

s-s' = o + (e 2 ) . (180) 

By substituting these expansions in Eqs. (|168|) and (|170|) . after straightforward calculations the asymptotic solutions 
(fTT!?]) and (TIT!))) follow identically. 
Q.E.D. 



We notice that the asymptotic expansion of the self-potential illustrated in THM.7 is required to reduce the non- 
local dependencies which are implicit in the definition of the fluid fields through the KDF. On the other hand, within 
the approximation obtained in THM.7 for the 4- vector K^, the RR equation (|165[) reduces to a local third-order 
ordinary differential equation. In particular, Eq. (ll70[) in THM.7 represents the analogue of the LAD equation for the 
single-particle dynamics, which contains the first derivative of the particle 4-acceleration (see also related discussion 
in Paper I). In view of this similarity, the asymptotic solution (|175l) can be further simplified adopting a second 
reduction-step of the same kind of that which leads to the LL form of the self- force for single charged particles [ll[ . 
This is obtained by assuming that the RR effect is only a small correction to the motion of the fluid. As a consequence, 
an iterative approximation can be adopted which permits to represent the self-force in terms of the instantaneous 
fluid forces. The latter include both the external EM field and the pressure forces. In particular, according to this 
method, to leading-order for the fluid 4-acceleration we have 

D = v/4 _ 1 

Ds (ext ^ " n ^ K ' 

where, for brevity we have introduced the notation 

F- t) = _JL_F^>- (i 82) 

The iteration gives 

+ -d ti P^U L d l Inn - -U l did,P^. (183) 
n n 

Substituting this expansion in Eq. ([175|) and invoking the symmetry property of the Faraday tensor provides for the 
first-order term the following approximation: 



K. 



' " + 3m„c 2 " + 3 



Mis — 2 . 

m c z l a 



m c z ^ n p 



t$> + M> , (184) 



where the first term on the rhs represents the mass-renormalization contribution, and hffl denotes the 4-vector 

h W = dl F { ;fu»U l - -A-F^F^Ut + S-s (T^U 1 ) (F {ext)kv U,) U„ (185) 

TYlriC TYlnC \ / \ ' 



while hp is given by 



= q ^-F ( *; t) d l P l ? + -d„P;U l d l \nn--U l d l d v Pr i + (186) 



+ ^-U li U k Ftfd l P lfi - -U, l U k d u P^U l d l Inn + -U„U k U l d t d v P^ , 
m c n H n n 

Eq. (jl84j) represents the fluid analogue of the LL approximation of the self- force holding for single particle dynamics, 
with the mass-renormalization term retained. In particular here we notice that: 

1) Eq. (|184p provides a local approximation of the fluid self- force carrying the contribution of the RR effect. In 
contrast to Eq. (|175p . thanks to the iterative reduction procedure only second-order derivatives of the position vector 
appear in this approximation. 

2) For consistency, Eq. (|184[) must be evaluated adopting the asymptotic expansion (jl T2[) also for the evaluation of 
the self-potential entering the definition of the fluid fields through the canonical momenta P„ in the KDF. 
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3) Moreover, consistent with the approximation in which the RR self-potential is small with respect to the external 
EM potential, also the asymptotic approximation (j!47p can be adopted, which allows one to treat explicitly in an 
asymptotic way all the implicit RR contributions. 

4) Finally, collecting together the analytical approximations provided by Eqs. (|147[) , (|172j) and (|184p . the fluid 
equations are reduced to a set of asymptotic local second-order PDEs. This provides a convenient representation also 
for Eulerian implementation schemes of the same equations. 

The detail comparison of Eqs. (|183l) - (|186l) with the literature is discussed in the next section. 



A. Retarded-time asymptotic expansion 

Despite the previous considerations, it is worth pointing out that, formally also for the fluid equations, an analogous 
result to THM.7 can be given. This is based on performing a Taylor expansion of the fluid RR force based on the 
retarded-time approximation. In this case, it is found that Eq. (|168p is approximated as 



-j(self) _ -j(self) 

A fi — A n 



where the retarded-time leading-order contribution A ' 

_( se Z/) 



{self) 
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[I + 0(e)], 
is simply given by 



while Eq. (|170p for the self- force, written in terms of defined in Eq. ([174p . becomes 

K»=K»\ t , [I + 0(e)], 
where the retarded-time leading-order contribution K^\ s , is now given by 
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with g' denoting here the 4-vector 
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(187) 



(188) 



(189) 
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(191) 



This alternative expansion has the distinctive advantage (with respect to the present-time expansion) of retaining 
all the physical properties of the exact fluid equations for the treatment of RR delay-time effects. This alternative 
formulation is relevant for comparisons with the point-particle treatment. 



XI. DISCUSSION AND COMPARISONS WITH LITERATURE 



In this section we analyze in detail the physical properties of the kinetic and fluid theory developed for the EM RR 
proble m, p roviding also a comparison with the literature. This concerns, in particular, the recent paper by Berezhiani 
ct al. 17], where an analogous research program is presented for the relativistic hydrodynamics with RR based on 
the LL solution of the self-force. 



A. Kinetic theory 



Let us start by considering the kinetic theory. The solution here obtained has the following key features: 

1) The kinetic theory adopts the Hamiltonian formulation of the RR problem here developed. The result is based 
on the exact analytical solution for the EM self-potential of finite-size charged particles, obtained in Paper I and 
Appendix A. 

2) The kinetic theory is here developed for systems of charged particles subject to an external mean-field EM 
interaction and the RR self-interaction produced by the same particles. Due to the non-local property of the RR 
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interaction, the formulation of kinetic theory is non-trivial. For this purpose, in contrast to previous literature, 
an axiomatic formulation of CSM is adopted. Its key element is the introduction of a suitable definition for the 
Lorentz-invariant probability-measure in the particle extended phase-space. As a consequence, the corresponding 
Liouville-Vlasov kinetic equation with the inclusion of the exact RR effect is achieved in Hamiltonian form, namely 
in such a way to preserve the phase-space canonical measure. For comparison, instead, previous literature approaches 
dealt with measure non-preserving phase-space dynamics. 

3) In particular, the kinetic theory has been developed within the canonical formalism representing the KDF in 
terms of the canonical state y = (r^,P^). For reference, in Appendix B the connection with the corresponding 
non-canonical treatment is provided. This in turn implies that non-local contributions associated to the self-potential 
(|i?2"j) enter implicitly in the definition of the corresponding fluid moments (|141[) - (I143[) . This is made possible only 
within the framework of the present exact formulation, in which the analytical solution for the self-potential is by 
construction non-divergent. This feature departs from recent approaches where instead non- Hamiltonian formulations 
were adopted, based on the LL point-like approximation of the RR self-force. In such a case in fact, the explicit 
dependence of the KDF in terms of the EM self-potential cannot be retained. 

4) Both the RR equation for single-particle dynamics and the kinetic equation for the KDF are of delay-type, 
reflecting the characteristic nature of the RR phenomenon. This property is completely missing from the previous 
literature on the subject, exclusively based on the LL local asymptotic approximation. 



B. Fluid theory 



For what concerns the fluid treatment, we notice that: 

1) Both the fluid fields and the fluid moment equations retain the standard form (available in the absence of RR 
effects) and can be equivalently represented in Eulerian or Lagrangian form. This follows from the exact representation 
here adopted both for the RR self-potential and the RR self-force. In both cases the only non-local dependencies are 
those associated to the position 4-vector. 

2) The exact fluid equations with the inclusion of the RR effect are delay-type PDEs. Because of this feature, their 
natural representation appears to be the Lagrangian form. In fact, the integration along the LPs must be in principle 
performed taking into account the retarded RR interaction. 

3) From the exact theory presented here it follows that each fluid equation of a given order does not depend 
on fluid fields of higher orders. For example, the momentum equation contains only second-order tensor fields, 
identified respectively with the plasma stress-energy tensor and the EM Faraday tensor. This result contrasts with 
the treatment given in Ref . (l7j where instead the asymptotic formulation based on the LL equation leads to moment 
equations involving higher-order tensor fields (for comparison, see also the related discussion in Appendix B). 

4) If a kinetic closure is chosen, then the fluid moments appearing in the fluid equations are all uniquely determined. 
In particular, the stress-energy tensor is prescribed in terms of the KDF. This implies that both implicit and explicit 
contributions of the RR effect appear in the resulting equations, carried respectively by the fluid fields and the EM self- 
force in the momentum equation. Remarkably, kinetic closure is achieved prescribing solely the pressure contribution 
carried by the stress-energy tensor. Instead, in the approach of Ref.[17| the closure conditions involve generally also 
the specification of higher-order moments of the KDF. 

5) An important feature of the exact fluid equations here obtained is that they can in principle be exactly imple- 
mented numerically adopting a Lagrangian scheme. 

6) A remarkable aspect of the present theory is that the relevant asymptotic expansions are performed only "a 
posteriori" after integration over the velocity space. This means that the approximations involved are introduced only 
on the configuration space- variables (i.e., the fluid fields) and not on the phase-space KDF. In particular, a convenient 
approximation is the one obtained in the short delay-time ordering, which reduces the non-local dependencies to local 
terms. As a consequence, the introduction of higher-order moments is ruled out by construction. 



C. Comparison with point-particle treatments 

The relevant comparison here is represented by Ref.JlT}. Such an approach is based on the adoption of the LL 
equation for the single-particle dynamics for the construction of the relativistic Vlasov-Maxwell description. The 
corresponding moment equations can be in principle adopted for the construction of a closed set of fluid equations. 
This requires however the specification of suitable closure-conditions. Let us briefly point out the novel features of the 
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current treatment for what concerns the adoption of the finite-size particle model in the construction of the kinetic 
and fluid descriptions. In detail: 

1) Both in the kinetic and fluid treatments the RR force is taken into account by means of a non-local interaction. 
This is an intrinsic feature of the assumed finite extension of the charged particle. In the fluid treatment, in particular, 
as shown above, the RR force can be parametrized in terms of the past Lagrangian fluid velocity and position. This 
permits to treat consistently the causal delay-time effects due to the finite-size of the particles. 

2) In validity of the asymptotic ordering given by Eq. (l99[) . an asymptotic retarded-time Hamiltonian approximation 
of the RR force based on a retarded-time expansion has been given for the fluid equations. This approximation 
preserves the basic physical features of the solution based on the exact form of the RR self-force. 

3) If the present-time asymptotic expansion is performed on the exact fluid moment equations, the resulting ex- 
pression of the fluid RR force obtained adopting the finite-size charge model appears different from that given in 
Ref.[I3. 

These conclusions enable us to carry out a detailed comparison with the literature, emphasizing the basic differences 
between kinetic and fluid treatments based on finite-size and point particles. 

A) Kinetic theory. 

The kinetic equation adopted in Ref.[ljJ is based on the LL equation (see therein Eqs.7 and 8). This means that 
the RR force in this approximation is non-conservative, non-variational and therefore non-Hamiltonian. In addition 
the LL equation: 1) does not retain finite delay-time effects characteristic of the RR phenomenon; 2) is not valid 
in the case of strong EM fields, where the iterative reduction scheme on which it is based, may fail; 3) ignores 
mass-renormalization effects (which are incompatible with the point-particle model). In contrast, the treatment of 
the relativistic Vlasov kinetic equation obtained here (see the Eulerian and Lagrangian equations (|I37[) and (|138[1 ) 
is qualitatively different. In fact, even if the resulting RR equation remains a second-order ODE, it is conservative, 
variational, Hamiltonian and applies for arbitrary external EM fields. Further remarkable aspects are related to the 
adoption of the finite-size charge model, in which the charge and mass distributions have the same support. As a 
consequence, in this case the self 4-potential is everywhere well-defined, contrary to the point particle model. In 
addition, this is prescribed analytically (see Appendix A), a feature which allows one to treat consistently the RR 
delay-time effects. 

B) Fluid theory. 

The fluid treatment here obtained is provided by the Eulerian Eqs. (|144[) - (|145|) or the equivalent Lagrangian equa- 
tions (|164j) - (ll65[) . The latter, considered as fluid equations, are manifestly not closed. However, the Hamiltonian 
formulation achieved here and holding for finite-size particles allows one to achieve a physically consistent kinetic 
closure condition, by prescribing uniquely the pressure tensor in Eq. (|165|) . We stress that in our treatment no 
higher-order moments need to be specified. In contrast, the corresponding Euler equation reported in Ref. [l7| (see 
Eqs.ll and 12 therein) actually depends also on a third-order tensor moment, which must be prescribed (see comments 
in Sec.IIIA of Ref. [17]). Let us now consider the asymptotic fluid treatments based on the present theory. These 
can be achieved invoking either the present-time or the retarded-time asymptotic expansions (see Section X). The 
first expansion is mostly relevant for comparisons with Ref.[l7j (given in THM.7) and enables one to achieve a local 
approximation of the delay-time effects carried by the RR force. However, remarkably, the resulting asymptotic fluid 
equations (|183l) - (|186|) remain qualitatively different from the corresponding ones given in Ref.pj}. In particular: 1) no 
higher-order moments appear after performing the Taylor expansion and the iteration scheme discussed after THM.7; 
2) a non-vanishing mass-correction contribution is now included (see first term on the rhs of Eq. (|184[ l). Finally, we 
mention that the retarded-time asymptotic expansion given by Eqs. (|187[) - (I191[) provides a novel approximation which 
retains basic properties of the exact solution. In particular: 1) it only applies for finite-size particles; 2) it relies on the 
Hamiltonian formulation of the RR problem and of the Vlasov-Maxwell treatment; 3) it permits to retain transient- 
time and delay-time effects; 4) it takes into account retarded mass-correction effects; 5) in this approximation the 
natural fluid description is Lagrangian. 



XII. CONCLUSIONS 

In this paper, novel results have been obtained concerning the kinetic and fluid descriptions of relativistic collisionless 
plasmas with the inclusion of EM RR effects. 

Relevant consequences of the variational form of the EM RR equation previously achieved for classical finite-size 
charged particles have been investigated. It has been shown that the non-local RR problem admits both Lagrangian 
and Hamiltonian representations in standard form, defined respectively in terms of effective Lagrangian and Hamil- 
tonian functions. A remarkable novel feature of the theory concerns the development of a Hamiltonian retarded-time 
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expansion of the RR force, which applies in validity of the short delay-time asymptotic ordering. On such a basis, 
the axiomatic formulation of classical statistical mechanics for relativistic collisionless plasmas with the inclusion of 
non-local RR effects has been presented. As a major result, the kinetic theory for such a system has been formulated 
in standard Hamiltonian form. The Liouville-Vlasov equation has been proved to hold in the extended phase-space, 
subject to non-local RR self-interactions. Remarkably, the non-local effects have been proved to enter the relativistic 
kinetic equation only through the retarded particle 4-position. As a consequence, the corresponding fluid moment 
equations can be determined in standard way by integration over the space of canonical momenta and cast both in 
Eulerian and Lagrangian forms. It has been pointed out that the exact relativistic fluid equations are intrinsically 
of delay- type and contain both implicit and explicit non-local contributions associated to the RR effect. The issue 
concerning the problem of fluid closure conditions has been discussed. In contrast with previous literature, it is found 
that in the present approach the closure conditions remain the standard ones, i.e., as in the absence of RR effects. 
Hence, the specification of higher-order moments of the KDF, for a given moment equation, is not required. Fi- 
nally, appropriate approximations have been obtained for the fluid equations by employing "a posteriori" the relevant 
asymptotic expansions applicable in the short delay-time ordering. This allows one to reduce the exact non-local 
equations either to a set of local PDEs or to retarded PDEs still retaining finite delay-time effects. 

The theory here developed has potential wide-ranging applications which concern the study of relativistic astrophys- 
ical plasmas for which RR emission processes are important. This involves, for example, plasmas in accretion disks, 
relativistic jets and active galactic nuclei. Other possible applications are also suggested for the case of laboratory 
plasmas generated in the presence of pulsed-laser sources. 



XIII. APPENDIX A: INTEGRAL REPRESENTATION FOR A 



{self) 



CASE OF A NON-ROTATING 



SPHERICAL-SHELL CHARGED PARTICLE 



In this Appendix we determine explicitly the integral representation of A^ el ^ for a non-rotating finite-size charged 
particle described by the model introduced in Paper I. We first remark that Eqs. (j6"2l and can be written for a 
spherically-symmetric charged particle of radius a > as 



^(s) = 0. 



(192) 
(193) 



Eq. (ll92j) defines the boundary dft on which the charge and mass are uniformly distributed, while Eq. (|193p represents 
the constraint of rigidity for the finite-size particle. We can use the information from Eq. (|192[) to define the internal 
and the external domains with respect to the mass and charge distributions. In particular, in terms of the generic 
displacement 4- vector X^ 6 M 4 defined as 



and subject to the constraint 



the following relations hold: 



X» Ufl (s) = 0, 

X^X^ < — er 2 '.external domain, 
X^X^ > —a 2 '.internal domain, 
X^X^ = = — er 2 '.boundary location. 



it is sufficient to know the solution of A 



(194) 



(195) 



(196) 



(self) 



As proved in Ref.[Tj|, for the evaluation of the action integral Sq 
the external domain. In this domain the EM self 4-potential generated by the non-rotating finite-size particle must 
necessarily coincide with that of a point particle carrying the same total mass and charge. The retarded 4-potential 
of a point charge represents a well-known result in the literature 
Green function approach. 



22] 



self-potential A 
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This can be easily obtained by means of the 
In particular, introducing the retarded Green function of a point particle G{r — r'), the 



takes the form 



A^n^r) = / d 4 r'G(r - r')f{r'), 



(197) 
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where G(r — r') is symmetric in both r and r' , is non-vanishing only for r° < r and satisfies in this domain the 
boundary-value problem 



UG{r - r') = 5^ (r - r') . 
G(r - r') = 0. 



As a consequence, written in integral form, the self-potential becomes 



A^ self ^(r,q) =2q d^<5(^ Q i? a ), 



with 



R° 



r a (s)- 



(198) 



(199) 



(200) 



This solution, derived for a point charged particle, also holds for the rigid finite-size non-rotating spherical shell in 
the external domain defined in Eq. (|196l) . As can be seen, this coincides with the result in Ref. , where a complete 
covariant solution for the EM self 4-potential A*jf el ^ holding in both internal and external domains has been obtained 
by adopting a derivation based on the principle of relativity and analogous to that outlined in Ref. fill ] for the point 
charge case. Notice that, contrary to the case of point particle, the self-potential (|199|) is well- defined also on the 
support of the charge, namely the ensemble on which the charge is distributed. 



XIV. APPENDIX B: NON-CANONICAL REPRESENTATION 



In this appendix we present the equivalent representation of the kinetic theory developed in section adopting non- 
canonical variables. For dcfinitcncss, let us introduce an arbitrary non-canonical phase-space diffcomorphism from T 
to r w , with r w denoting a transformed phase-space having the same dimension of T, 



y = (r' I ,P M )^w = w(y), 



(201) 



where, for example, w can be identified with the non-canonical state y nc = (r^,p^) defined in Eq. (|147|) or with 
y u = (r / - l ,u fi ). In the second case the transformation, following from Eq. (|96|) . is realized by 



Up = Pfj, - - 
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The transformed RR equation in the variables y u becomes therefore: 



ds 
du^ 

ds 
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where F u = if^u 

r or 



Denoting now by 



A(w(s)) = 



dy(s) 
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(204) 
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the KDF mapped onto the transformed phase-space T w by the KDF / (y (s)), the differential Liouville-Vlasov equation 
(|127p requires 



d 
ds 



<9w (,s) 
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/i(w(s)) 



0. 



where Wo = w (so). At the same time, Eq. (|127l) also implies, thanks to the chain rule: 

d 



ds 



/(y(w(,s))) = 0, 
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which for consistency delivers the well-known differential identity 



Prom Eq. pTJTj) it follows 
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This equation can be represented, for example, in terms of w = y u . In this case, due to the chain rule 
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while, thanks to Liouville theorem 
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As an application of the result, it follows that, if the LL approximation is introduced for the 4-vector F^, namely 
Eqs. (|204|) and (|205|) are replaced with asymptotic equations of the form 



l LL> 



dr>* 
LL 

ds 

du LL _ P M 

ds ~ LL ' 



(213) 
(214) 



where F£ L is the total EM force in this approximation, then Eq. (|210p recovers the expression reported in Ref. [17|. 
This provides the connection with the exact canonical theory here developed. We remark, however, that since the 
LL equation is only asymptotic, the mapping between the canonical state y = (r M , P^) and yll = ( r LL^ u ll^l) is also 
intrinsically asymptotic. Therefore, Eqs. (|213[) and (|214p remain necessarily non- variational and non-canonical. 
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